First Study for the Pentaquark Potential in SU(3) Lattice QCD 
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The static penta-quark (5Q) potential Vbq is studied in SU(3) lattice QCD with 16 3 x 32 and 
/3=6.0 at the quenched level. From the 5Q Wilson loop, Vsq is calculated in a gauge- invariant 
manner, with the smearing method to enhance the ground-state component. Vsq is well described 
by the OGE plus multi-Y Ansatz: a sum of the OGE Coulomb term and the multi-Y-type linear 
term proportional to the minimal total length of the flux-tube linking the five quarks. Comparing 
with QQ and 3Q potentials, we find a universality of the string tension, <7qq ~ 03Q ~ 05Q, and the 
OGE result for Coulomb coefficients. 
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The inter-quark force is one of the elementary quan- 
tities for the study of the multi-quark system in the 
quark model. As for baryons, our group recently studied 
the three-quark (3Q) potential V3Q in detail with lattice 
QCD, and clarified that it obeys the Coulomb plus Y- 
type linear potential^]. However, no one knows the inter- 
quark force from QCD in the exotic multi-quark sys- 
tem such as tetra-quark mesons (QQ-QQ), penta-quark 
baryons (4Q-Q), dibaryons (6Q) and so on. 

Very recently, an exotic anti-strange baryon O + (1540) 
with S = +1 was experimentally discovered at SPring- 
8 (LEPS), and was confirmed by ITEP (DIANA), 
JLab(CLAS) and ELSA(SAPHIR)0. The 6+(1540) 
was theoretically predicted in the Skyrme model[3j, and 
is regarded as a penta-quark (5Q) baryon of u 2 d 2 s in the 
valence-quark picture. Another 5Q baryon S~~(1862) 
was found by CERN(NA49)0, and also an anti-charmed 
5Q baryon 9 C (3099) was found by HERA(H1)0. 

Accordingly, many theoretical studies^] have been 
done for the 5Q baryon using various approaches such 
as lattice QCD0J3, the constituent quark modelQ, the 
diquark model[l(|, the QCD sum rule ( llj . the flux-tube 
model0, the string theory and so on. However, 
there are several puzzling problems on the 8 + (1540): its 
mass seems to be rather small and its decay width is ex- 
tremely narrowj2J. To solve them, one encounters the 
many-body problem of quarks, and therefore it is quite 
desired to clarify the inter-quark force in the multi-quark 
system based on QCD. 

In this paper, motivated by the recent discovery of 
the penta-quark baryons, we perform the first study of 
the static penta-quark (5Q) potential V5Q, i.e., the inter- 
quark force in the 5Q system, in SU(3) lattice QCD with 
(3=6.0 and 16 3 x 32 at the quenched level. Note that the 
lattice QCD result of V5Q presents a key information in 
modeling the multi-quark system based on QCD. 

For the penta-quark system, we investigate the QQ-Q- 



QQ type configuration with the two "QQ clusters" be- 
longing to the 3* representation of the SU(3) color as 
shown in Fig.l, since this type of the 5Q configuration is 
expected to have a small energy and seems to be natural 
as a realistic candidate of the 8 + (1540). Indeed, in the 
perturbative sense, an attractive (repulsive) force acts be- 
tween two quarks, when their total SU(3) color belongs 
to 3* (6). Therefore, the nearest QQ cluster tends to 
form 3* rather than 6 in the low- lying 5Q system, which 
leads to the 3*-diquark model|lO|. 




FIG. 1: The QQ-Q-QQ type configuration for the penta- 
quark system. The two QQ clusters belong to the 3* repre- 
sentation of the color SU(3). 



Similar to the derivation of the Q-Q (3Q) potential 
from the (3Q) Wilson loop, the 5Q static potential Vj>q 
can be calculated with the 5Q Wilson loop W5Q, which 
is defined in a gauge- invariant manner as shown in Fig. 2. 

We define the 5Q Wilson loop W 5Q as 

W 5Q = i e Qbc e a ' fc ' c 'M aa '(L 3 L 12 L 4 ) bb '(i? 3 i?i 2 i?4) cc ', (1) 
where M, L i: R\(i = 1, 2, 3, 4) are given by 



M, U, Ri = Pexp{i«? / dx^A^x)} € SU(3) C . (2) 

JM,Li,Ri 

As shown in Fig. 2, M, Li, Ri(i = 3,4) are line-like vari- 
ables and Li,Ri(i — 1,2) are staple-like variables. Here, 
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FIG. 2: The penta-quark (5Q) Wilson loop W 5 q for the 5Q 
potential Vsq. The contours M, Li, Ri(i = 3,4) are line-like 
and Li,Ri(i — 1,2) are staple-like. The 5Q gauge-invariant 
state is generated at t = and is annihilated ait — T with the 
five quarks (4Q-Q) being spatially fixed in R 3 for < t < T. 



Li2,Ru G SU(3) C arc defined as 



t a' a 1 abc a'b'c 

Ll2 = 2 6 



/ bb' J cc 



T>(1 CI 

H 12 



_^abc^a b i 

2 
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.(3) 



Note that the 5Q Wilson loop Wsq is gauge invariant, 
and its gauge invariance is owing to the nontrivial as- 
signment of the color indices of L 12 and R\ 2 in Eq.Q. 
(Recall that the "two quark lines" combining into the 3* 
representation correspond to an "antiquark line" as the 
color current.) 

In principle, the ground-state 5Q potential Vsq is 
obtained from the 5Q Wilson loop (W^q) as V5Q = 
— limT^oo 4; ^(Wsq). However, the practical lattice cal- 
culation is done with a finite region of T, where excited- 
state contributions remain. For the accurate measure- 
ment of V5Q in lattice QCD, we use the gauge-covariant 
smearing method 0] to enhance the ground-state compo- 
nent of the 5Q state in the 5Q Wilson loop. 

The smearing is known to be a powerful method 
for the accurate measurement of the Q-Q and the 3Q 
potentials^], and is expressed as the iterative replace- 
ment of the spatial link variables Ui(s) (i=l,2,3) by the 
obscured link variables Ui(s) G SU(3) C which maximizes 
Retr {Uj(s)Vi(s)} with 

Vi(s) = aU^s) + ]T ]T{[/±, (*){/,(* ± (s + I)} (4) 



We 



with the simplified notation of U-j = Uj(s — j). 
here adopt a — 2.3 and the iteration number iV smr = 40, 
which lead to a large enhancement of the ground-state 
component in the 5Q Wilson loop at /3=6.0. 

Now, we proceed the actual lattice QCD calculation 
for the 5Q potential V5Q 0. We generate 150 gauge con- 
figurations using SU(3) C lattice QCD with the standard 
action with (3 = 6.0 and 16 3 x 32 at the quenched level. 
The gauge configurations are taken every 500 sweeps af- 
ter a thermalization of 5000 sweeps using the pseudo- 
heat-bath algorithm. The lattice spacing a is estimated 
as a ~ 0.104fm from the string tension <r=0.89 GeV/fm 
in the Q-Q potential Vqq 0. 



As for the 5Q configuration, we consider the QQ-Q- 
QQ type configuration as shown in Fig.l. Note that the 
multi-quark system including four or more quarks can 
take a three-dimensional shape, while the QQ and the 3Q 
systems can take only planar configuration!^, ^2- Then, 
we investigate both the planar 5Q configuration as shown 
in Fig. 3 and the twisted 5Q configuration as shown in 
Fig. 4. In this paper, we take d\ = d 2 = (I3 = d± = d, 
and present the lattice QCD result of V5Q in terms of 
(d, hx,h 2 ). 
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FIG. 3: A planar configuration of the penta-quark system. 
Q1Q2 is parallel to Q3Q4, and H1H2 is perpendicular to Q1Q2 
and Q3Q4. Here, we take di — c?2 = d$ — di = d. 

Q, 




FIG. 4: A twisted configuration of the penta-quark system. 
Q1Q2 is perpendicular to Q3Q4, and H1H2 is perpendicular 
to Q1Q2 and Q3Q4. Here, we take di — c?2 = dz — di = d. 



For these types of 5Q configurations, we calculate 
the 5Q potential V5Q from the 5Q Wilson loop (Wsq) 
using the smearing method. Owing to the smearing, 
the ground-state component is largely enhanced, and 
therefore the 5Q Wilson loop (Wsq) composed with the 
smeared link variable exhibits a single-exponential be- 
havior as (W5Q) ~ e~ V5C>T even for a small value of T. 
Then, for each 5Q configuration, we extract Vsq from the 
least squares fit with the single-exponential form 



(W 5 Q> = Ce~ v ^ T 



(5) 



in the range of T u 
prefactor C physically means the ground-state overlap, 
and C ~ 1 corresponds to the quasi-ground-state. Here, 
we choose the fit range of T such that the stability of the 
"effective mass" V(T) = ln{(W 5Q (T)) / (W 5Q (T + 1))} is 
observed in the range of T min < T < T max — 1 . For the 
lattice calculation of (Wsq), we use the translational and 
the rotational symmetries on lattices. 

For 56 different patterns of the 5Q configurations as 
shown in Figs. 3 and 4, we present the lattice QCD data 
for the 5Q potential Vj>q together with the ground-state 
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TABLE I: Lattice QCD results for the penta-quark potential 
V5Q for the planar 5Q configuration labeled by (d,hi,h2) as 
shown in Fig. 3. We list also the ground-state overlap C, the 
fit range of T and the theoretical form Vgn 00 * of the OGE plus 
multi-Y Ansatz 10 with (Asq^sq) fixed to be (A3q,<j3q) in 
V 3 q in Ref.Q- All the data are measured in the lattice unit. 



(d, hi,h 2 ) 


V 5Q 


C 


-Mnin Tin ax 


t /thcor 


(1,1,1) 


f.4452(ff) 


0.9539(21) 


2-7 


1.4433 


(1,1,2) 


f.5409(13) 


0.9506(25) 


2-8 


1.5414 


(1,1,3) 


L6f77(f9) 


0.9512(33) 


2-7 


1.6146 


(1,1,4) 


1.6793(20) 


0.9431(35) 


2-7 


1.6767 


(1,1,5) 


1.7381(23) 


0.9394(40) 


2-6 


1.7332 


(1,1,6) 


1.7918(28) 


0.9311(49) 


2-6 


1.7866 


(1,1,7) 


1.8441(31) 


0.9232(56) 


2-6 


1.8380 


(1,2,2) 


1.6314(17) 


0.9503(29) 


2-6 


1.6322 


(1,2,3) 


1.7011(20) 


0.9427(34) 


2-5 


1.7021 


(1,2,4) 


1.7680(25) 


0.9478(43) 


2-6 


1.7623 


(1,2,5) 


1.8190(29) 


0.9297(50) 


2-6 


1.8177 


(1,2,6) 


1.8717(33) 


0.9205(57) 


2-5 


1.8704 


(1,3,3) 


1.7723(24) 


0.9405(39) 


2-4 


1.7702 


(1,3,4) 


1.8336(28) 


0.9351(49) 


2-7 


1.8293 


(1,3,5) 


1.8913(33) 


0.9320(62) 


2-5 


1.8839 


(1,4,4) 


1.8939(31) 


0.9293(53) 


2-4 


1.8877 


(2,1,1) 


1.7531(23) 


0.9393(40) 


2-5 


1.7515 


(2,2,2) 


1.8803(31) 


0.9292(54) 


2-6 


1.8887 


(2,3,3) 


2.0030(37) 


0.9284(64) 


2-5 


2.0098 


(2,4,4) 


2.1122(49) 


0.9116(91) 


2-5 


2.1211 


(3,1,1) 


1.9734(37) 


0.9138(66) 


2-5 


1.9850 


(3,2,2) 


2.0811(45) 


0.9070(75) 


2-5 


2.0942 


(3,3,3) 


2.1886(53) 


0.9003(92) 


2-4 


2.2047 


(3,4,4) 


2.3043(68) 


0.9084(113) 


2-5 


2.3105 


(4,1,1) 


2.1697(60) 


0.8948(100) 


2-5 


2.1958 


(4,2,2) 


2.2734(60) 


0.8890(100) 


2-5 


2.2829 


(4,3,3) 


2.3657(73) 


0.8606(120) 


2-4 


2.3864 


(4,4,4) 


2.4706(104) 


0.8534(164) 


2-5 


2.4884 



overlap C in Table I and II. The statistical errors are 
estimated with the jackknife method. We find a large 
ground-state overlap as C > 0.85 for almost all 5Q con- 
figurations. 

Next, we consider the theoretical form of the 5Q po- 
tential V5Q. The lattice QCD studies^] at the quenched 
level show that the Q-Q potential Vqq takes a form of 



V Q d(r) = - 



■ + °w 

and the 3Q potential V3Q takes a form of 

1 



Cqq, 



V f 



3Q 



a 



3Q, 



(0) 



(7) 



where L m ; n denotes the minimal value of total length of 
color flux tubes linking the three quarks. In fact, both 



TABLE II: Lattice QCD results for the penta-quark potential 
V5Q for the twisted 5Q configuration labeled by (d, hi, /12) as 
shown in Fig. 4. The notations are the same in Table I. 



\a, rii , ri2 ) 






T ■ T 

-*- mm J- max 


r/thcor 
^5Q 


( 1 1 1 A 


1 a A7«f r >'i\ 
1.44 (0{Z4 ) 


fl QQ7G/' fi/l ^ 

U.yo / 0(^04 J 


O-O 


1 a a p;fi 
i.440o 


h 1 


1 X>A QQ/'l A \ 
1.04oo^l4j 


u.yozo^zo ) 


O fi 
Z-o 


1 KA 1 O 

i.04iy 


fi 1 

IM>°J 


X.OiOO^l / ) 


n Q/i p»Qf qi ^ 
u.y4oy^oi ) 


O 

Z-0 


1 fil A fi 
1.014O 


( 1 1 A \ 

IMA) 


1 fi7fi7 1 01 \ 
l.D/D/ {11 ) 


fl OQ7A(' Q7^ 

u.yo (\j{4 ( ) 


O fi 
Z-0 


1 fi7fi7 
1.0 / O / 


( 1 1 ^\ 
IMA) 


1 7QfiP;/'00 s \ 
1. ( 400{ZZ ) 


Q1Xi7( AO\ 

u.yoo / ^4z ) 


^ 
Z-0 


1 7QQ9 
1. / OOZ 


IMA) 


1 7Gio/'ofi s \ 
1. f yiz^zo ) 


fl QOQ7^ AO\ 

u.yzy / ^4z ) 


1 

Z-4 


1 7fififi 
1. / oOO 




1 fiQQ7/'fifi s \ 
l.OOO I {OO ) 


U.oyoo\Zol ) 


o-O 


1 fiQfifl 
I.000U 


{1,2,2) 


1 fi^no { 1 f<\ 

l.OoUZ^lO ) 


Q/170/'0Q > \ 

u.y4 / z^zy ) 


9 fi 
Z-o 


1 fiQ9/l 
1.0oZ4 


{1,2,4) 


1 Tnoo/'i q\ 

1. / UZZ^lo ) 


n Q/i a Kt qo^ 

U.y440^oZ J 


1 

Z-4 


1 7099 
1. / UZZ 


IM) 4 J 


x. ( OO / ^ZO ) 


u.y^z / 1 i ±'-tL 1 


Z-0 


1 7fi9/1 
1. ( OZ4 


{l,2,o) 


l.oZoZ^OU ) 


u.yooo^oi ) 


O 7 
Z- ( 


1 fil 77 
l.ol / / 


{l,2,b) 


1 fi70fi/'Q0\ 
l.o / Zo{42 ) 


u.yzou^oo ) 


O fi 
Z-o 


1 fi7fl/1 

l.o / U4 




1 771fV0/1 s \ 
1. 11 1U^Z4 ) 


u.yo / o^4z ) 


O fi 
Z-0 


1 7709 
1. / / UZ 


( 1 Q A \ 


1 fiQ0fi/07 s \ 
l.ooZO^Z ( ) 


u.yooo^o ) 


9 1 

Z-4 


1 Q90Q 

l.oZyo 


{1,4,0) 


1.0i)DZ{4Z ) 


u.yoy4^oo j 


O 7 
Z- / 


1 fifi^Q 

i.oooy 


(~\ A A} 

UA 4 J 


l.oyOU^oU ) 


n qqi 

u.yoio^oz ) 


O /I 
Z-4 


1 fifi77 
l.OO / ( 


\2,i,i) 


1 77Q^/'OQ\ 
l.(( oO^ZO ) 


u.yo / / ^41 ) 


O fi 
Z-0 


1 7fi1 K 

1. / 010 


( o\ 
{2,2,2) 


1 OQQO/Qn\ 

l.oooz(oU ) 


u.yz / y^o4 j 


O 7 
Z- / 


1 Qonn 
l-ooyy 


{2,4,4) 


nm 1 /'Q7 S \ 

Z.UU1 1 ^0 / j 


u.yzoo^04 j 


9 '■, 
Z-0 


ni nn 

Z.U1UU 




Oil ^^('/IQA 
z. i ioo^4y ) 


u.y 10 / ^oo ) 


9 fi 
Z-0 


01010 

Z. 1Z1Z 


(3,1,1) 


z.UU4y(o / ) 


U.yUoz^o / J 


Z-0 


nnno 
z.UUUo 


(3,2,2) 


2.0873(38) 


0.8987(69) 


2-4 


2.0973 


(3,3,3) 


2.1870(54) 


0.8912(92) 


2-6 


2.2055 


(3,4,4) 


2.3021(68) 


0.9019(113) 


2-5 


2.3108 


(4,1,1) 


2.2141(65) 


0.8741(107) 


2-6 


2.2155 


(4,2,2) 


2.2874(64) 


0.8768(107) 


2-5 


2.2879 


(4,3,3) 


2.3715(70) 


0.8577(118) 


2-4 


2.3880 


(4,4,4) 


2.4680(94) 


0.8459(149) 


2-4 


2.4890 



I^qq and V3Q are described by a sum of the short-distance 
one-gluon-exchangefOGE) result and the long-distance 
flux-tube result 0,131 . 

For the static penta-quark (5Q) system, we find that 
the lattice QCD results are well described by the OGE 
plus multi-Y Ansatz: a sum of the OGE Coulomb term 
and the multi-Y type linear term 



5Q 



9 2 V- 



i j 



4tt ^ \n 

i<3 



C 5 



with 



-Aq{( 
l , 1 



= r t 



ri2 
1 

ru 



1.1.1 1 

— ) + 7^(— + — + 

^34 2 r 15 r 25 

1 1 >, 

h + )} + <T5QLn 

f23 7*24 



1 



+ — ) 



1 

ill + C 5Q (8) 



I and ith quark location in Fig.l. 
Here, L m \ n is the minimal length of the flux-tube linking 
five quarks as shown in Fig.l. (For the extreme case, e.g., 
d > y/3hi, we here assume that the flux-tube is formed 
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as the two straight lines on Q1Q5 and Q2Q5, considering 
the color combination, although there may appear several 
possibilities as the "flip- flop" .) 

Note that there appear three kinds of Coulomb coef- 
ficients (A5Q, 5^5Q, j^4.5q) in the penta-quark system, 
while only one Coulomb coefficient, ^4qq or A 3 q, appears 
in the QQ or the 3Q system. Here, the Coulomb coeffi- 
cient A5Q in Eq.® corresponds to A 3 q or in terms 
of the OGE result. 

We add in Table I and II the theoretical form V^ eor of 
the OGE plus multi-Y Ansatz © with (A 5 q,(75q) fixed to 
be (^43q,o"3q) in the 3Q potential V3Q obtained in Ref.0], 
i.e., A 5Q = A 3 q ~ 0.1366, <t 5 q = a 3Q ~ 0.046a -2 and 
C5Q ~ 1.58a _1 . (Note that there is no adjustable param- 
eter for the theoretical form V 5 Q Cor besides an irrelevant 
constant C5Q, since A$q and CT5Q are fixed to be A 3 q 
and (T3Q, respectively.) Thus, the 5Q potential V5Q is 
found to be well described by the OGE Coulomb plus 
multi-Y-type linear potential. 

We show in Fig. 5 typical examples of the lattice QCD 
data for the penta-quark potential V5Q. The symbols de- 
note the lattice data, and the curves denote the theoreti- 
cal form of the OGE plus multi-Y Ansatz with (A5q,<75q) 
fixed to be (A 3 q,<j 3 q). One finds a good agreement be- 
tween the lattice QCD data and the theoretical curves. 
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FIG. 5: Lattice QCD results of the penta-quark potential 
V5Q for the planar 5Q configuration with hi = hi = h in 
Fig. 3 in the lattice unit. Each 5Q system is labeled by d and 
h. The symbols denote the lattice data, and the curves the 
theoretical form of the OGE plus multi-Y Ansatz. 

Note that the planar and the twisted 5Q configura- 
tions with the same (d, hi, h?) are almost degenerate, al- 
though the energy of the planar one is slightly smaller. 
In terms of the OGE plus multi-Y Ansatz, the only en- 
ergy difference between the two states originates from 
a small difference of the Coulomb interaction between 
Qi(i — 1,2) and Qj(j — 3,4), where the Coulomb coef- 
ficient is reduced as | j 4 5 q(~ ^Aqq). Then, no special 
configuration is favored in the 5Q system in terms of the 
energy. This fact also indicates that the 5Q system is un- 
stable against the twisted motion of the two QQ clusters 



as shown in Fig. 4. In fact, general 50 s yste ms tend to 
take a three-dimensional configuration |8|, ll2j m terms of 
the entropy. 

From the comparison with the QQ and the 3Q poten- 
tials, the universality of the string tension and the OGE 
result are found among QQ, 3Q and 5Q systems as 

CT QQ — CT 3Q — c 5 q, -Aqq~A 3 q~A 5 q. (9) 

This result supports the flux-tube picture on the confine- 
ment mechanism even for the multi-quark system 

To conclude, we have performed the first study of the 
penta-quark potential in lattice QCD, and have found 
that the 5Q potential is well reproduced by the OGE 
Coulomb plus multi-Y-type linear potential. 
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